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We consider the electronic properties of layered molecular crystals of the type 0-D2A where A is 
an anion and D is a donor molecule such as BEDT-TTF [where BEDT-TTF is bis-(ethylenedithia- 
00 ' tetrathiafulvalene)] which is arranged in the 9 type pattern within the layers. We argue that the 

simplest strongly correlated electron model that can describe the rich phase diagram of these ma- 
terials is the extended Hubbard model on the square lattice at a quarter filling. In the limit where 
the Coulomb repulsion on a single site is large, the nearest-neighbour Coulomb repulsion V plays a 
$H ■ crucial role. When V is much larger than the intermolecular hopping integral t the ground state is an 

' insulator with charge ordering. In this phase antiferromagnetism arises due to a novel fourth-order 

superexchange process around a plaquette on the square lattice. We argue that the charge ordered 
■ phase is destroyed below a critical non-zero value V , of the order of t. Slave boson theory is used 

to explicitly demonstrate this for the SU(N) generalisation of the model, in the large N limit. We 
also discuss the relevance of the model to the all-organic family /3"-(BEDT-TTF)2SF5YS03 where 
^ • Y = CH2CF2, CH2, CHF. 

o 

PACS numbers: 71.27.+a, 71.1Q.Fd, 74.70.Kn, 71.45. Lr 
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Layered organic |incjlecular crystals based on the BEDT-TTF molecule [where BEDT-TTF is bis-(ethylenedithia- 
tetrathiafulvalene)]Eria are rich examples of strongly correlated electron systems in low dimensions. The k-(BEDT- 
■ TTF)2X family has many similarities to the cuprate suuecconductors, such as the proximity of superconductivity 
"Y^ ' to an antiferromagnetic Mott-Hubbard insulating phase. It has recently been argued that the simplest strongly 
d ■ correlated |electron model that can describe this family is a Hubbard moded.x)n an anisotropic triangular lattice at 
^ . half-fiUing.El This model should also describe the /3-(BEDT-TTF)2X family.Lftl As the anion X or pressure is varied 
I ' the 6*- (BEDT-TTF) »X_family exhibits metallic, superconducting, insulating, antiferromagnetic, charge ordeiied, and 
spin gapped phasesE'EHl The family 6'-(BETS)2X [where BETS is bis-(ethylencdithio-tetraselenafulvalfine)]EJ is also 
^ • found to exhibit a metal-insulator transition with a transition temperature which varies with the anion.E3 The recently 
Q ; synthesized family /3"-(BEDT-TTF)2SF5YS03 where Y = CH2CF2, CH2, CHF has attracted considerable interest 
' from chemists because the anion is purely organic and Y = CH2CF2 is the first rawely organic superconductor.Ej 
Insulating, charge ordered and spin gapped phases are observed for different anionsJi-l 

Given the complexity of the details of the chemistry, crystal structure and band structures of these matecials it is 
^ important to define the simplest possible many-body Hamiltonian that can capture the essential physics .cJ This is 
similar in spirit to the way one studies the Hubbapd|raiid t — J models on a square lattice in order to understand the 
cuprate superconductoijSu. Several previous studiestSila of the metal- insulator transition and the magnetic properties 
of the insulating phasecZI of the 6 family have been interpreted in terms of the Hubbard model. Jlowever, this is 
inadequate because, at quarter filling, the Hubbard model is expected to be metallic. Following Seoll3, we argue that 
the nearest neighb our . Co ulomb repulsion V plays a crucial role in these materials. This has been emphasized before 
for other organics.BBEl'y 

Specifically, the simplest possible strongly correlated electron model that can describe the competition between all 
of the above phases is the extended Hubbard model (or a t — model with no double occupancy) at quarter-filling 
on a square lattice. We show that in the charge-ordered insulating phase antiferromagnetic interactions arise due to a 
novel fourth-order superexchange around a plaquette on the square lattice. We then consider the SU(N) generalization 
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of the t — V model and perform a slave boson study which becomes exact in the limit of large N. We find there is 
a critical value of the ratio V/t, above which the metallic Fermi liquid phase undergoes an instability to a charge 
ordered state. 



A. Review of experimental properties of 6'-(BEDT-TTF)2X 

Properties of the 6'-(BEDT-TTF)2X family have recently been reviewed by Mori et aL0@ and Seo[l. The arrange- 
ment of the BEDT-TTF or BETS molecules within a layer of the 6 structure is shown in Fig. |l|. For most anions, X, 
the materials undergo a metal-insulator transition at a temperature that decreases-svith increasing band-width; the 
latter is directly correlated with the angle between the molecules within the layerst^ (see Kg. ^). The temperature 
at which the metal-insulator transition occurs generally increases with increasiag pressure.li3 

6'-(BDT-TTP)2Cu(NCS)2 undergoes a metal-insulator transition at 250 K.l^ At low temperatures the charge gap 
deduced from the activation energy of the conductivity is about 100 meV. Raman-active vibrational modes (associated 
with the stretching of carbon-carbon double bonds) of the BPT-TTP molecules are sensitive to the charge on the 
molecule. In the insulating phase these modes split, consistent with charge ordering.Ej Evidence for short-range 
charge ordering along the c-axis dipeiction (the vertical direction in Fig. |l|) was found in the insulating phase of 
X=CsCo(SCN)4 by x-ray scatteringEj For X=CsZn(SCN)4, the principal axes of the g-tensor associated with electron 
spin resonance undergo a rotation at 20 K; this has been interpreted as a change in the electronic state^ For 
X=RbZn(SCN)4 there is a metal-insulator transition at 190 K; there is then a dimerization in the c-direction.li3 The 
magnetic susceptibility shows no features at this transition and between 50 and 190 K has been fit to that for a two- 
dimensional antiferromagnetic Heisenberg model with exchange J — 100 K. Below 50 K the susceptibility is consistent 
with a spin gap of about 4&iK. There is evidence from nuclear magnetic resonance for charge ordering below 190 K and 
of a spin gap below 10 For X= CuW|CN)[N(CN)2]2 there is a metal-insulator transition at 220 K; the charge gap 
at low temperatures is about 200 meV.EZI The magnetic susceptibility shows no features at this transition and between 
33 and 220 K has been fit to that for a two-dimensional antiferromagnetic Heisenberg model with J = 48 K. There 
is no sign of Neel order but below 30 K the susceptibility decreases rapidly, suggesting formation of a spin gap. The 
only member of the 6'-(BEDT-TTF)2X family that is superconducting is X=l3 which has a transition temperature 
of 3.6 K. The Fermi surface of the metallic phase has been mapped out using angular-dependent magnetoresistance 
and magnetic oscillations. Ej Several of the 0-(BETS)2X family undergo a metal-insulator transition and several are 
metallic down to 4 K (see Table |). 

6'-(CiTET-TTF)2Briwhere CiTET-TTF is bis(methylthio)ethylenedithio-tetrathiofulvalene] is an insulator with a 
charge gap of 600 meV.EZl The magnetic susceptibility between 8 and 290 K has been fit to that for a two-dimensional 
antiferromagnetic Heisenberg model with J = 6 K. Below 3 K the susceptibility depends on the field direction, 
suggesting the formation (j>fc|Neel order. 

As emphasized by Moric3, and illustrated schematically in Figure |2[ Table || shows the general trend that as the 
band width (which is roughly proportional to tp) increases the transition of the metal-insulator transition decreases. 

II. THE EXTENDED HUBBARD MODEL 

The arrangement of the BEDT-TTF molecules within a layer of 6'-(BEDT-TTF)2X is shown schematically in Figure 
1. Values of the intermolecular hopping integrals, calculated using the Hiickel approximation are given in Table |. If 
there is complete charge transfer of one electron onto each anion X there is an average of half a hole per molecule and 
so the electronic bands will beuauaptetJGilled. 

Band structure calculationsoEZlO'ESE^ predict that all these materials are metallic. Hence, the different phases 
must be due to strong electronic interactions. The Hubbard interaction U describing the Coulomb repulsion between 
two electroaSijOri the same BEDT-TTF molecule has been estimated from quantum chemistry calculations to be 
about 4 eV.EilE3 Thus U is much larger than the band width associated with the hopping integrals given in Table 1, 
confirming that these are strongly correlated materials. The Hubbard model on the square lattice at quarter filling is 
expected to be always metallic and so one must consider longer range Coulomb interactions to explain the existence of 
insulating, charge ordered and antiferromagnetic phases. The nearest neighbour Coulomb repulsion between yajioua 
arrangements of pairs of BEDT-TTF molecules has also been estimated from quantum chemistry calculatioufSOS 
and is generally found to have a value of about 2-3 eV. It is approximately given by Coulomb's law y ~ 14 eV/R, 
where R is in A. Mori calculated y as a function of the angle 9; variations of about ten per cent occur in the range 
(98 - 130 degrees) relevant to the 6'-type materials. These calculations involve isolated pairs of molecules and so one 
expects that the values of U and 1/ in a molecular crystal to be less than this due to screening. Hubbard has discussed 
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this for the case of TTF-TCNQ arguing that U and V are both decreased by a factor of about two.B Actually in 
Section [II, from experimentally determined charge gaps, we estimate values of V of the order of a few hundred meV. 

35^ 



In materials consisting of dimers of BEDT-TTF molecules the difference U — V can be estimated from the chara 
transfer excitation seen in optical absorption measurements. For (BEDT-TTF)IIgBr3 it is estimated to be 0.7 eV.I 
Thus we are led to the extended Hubbard model on the anisotropic triangular lattice at half filling. 

Table | shows that for many of the 6 materials, tc <SC tp and so, as a first approximation, we neglect the diagonal 
hopping tc- This means we are left with a square lattice model. In Section |^ we will show that this diagonal term 
has little effect on the metal-insulator transition. The Hamiltonian is 

H = t {c\^Cj„ + c]^c^a) + C/^ n,tn,x + y ^ rurij ~ ^ij^^^'^ (1) 

where U is the Coulomb repulsion between two electrons on the same site, V is the nearest neighbour Coulomb 
repulsion, < ij > denotes nearest neighbours, and /x is the chemical potential. 

If we consider the large U limit and so preclude doubly occupied sites, the Hamiltonian then reduces to 



where P projects out the doubly occupied states. We refer to this as the t — V model. 



For large V/t the ground state is an insulator with charge ordering (Section [H). We expect that for small V/t a 
metallic phase exists because the quarter- filled Fermi surface is poorly nested. Hence, as V/t decreases the charge 
ordering should be destroyed at a non-zero value of, V/t. We are not aware of any previous study of this model on 
the square lattice. Mazumdar, Clay, and CampbellEl have studied coupled chains of the extended Hubbard model at 
quarter filling. Most of their numerical results are for U = 6t|| and V|| = t|| (where the || refers to the chain direction) 
so that V^i is smaller than the critical value necessary to form the charge ordered state considered here. They find 
that when interactions with phonons are included there is a tendency to formation of a bond-otdcr wave. Henley and 
Zhang recently considered a similar model involving spinless fermions on the square latticeBj Mila considered the 
extended Hubbard model on the square lattice at a quarter filling with finite U and infinite VS3 In Section IV we 
consider the SU(N) generalization of the t — V model and show that in the large N limit, slave boson theory implies 
that there is a critical value of V/t above which the metallic phase becomes unstable to formation of charge ordering. 

We now briefly review previous work on the extended Hubbard model (in the large U limit) at quarter filling on 



different lattices. Numerical calculations s 
a charge-density wave insulator at = 2t. 



iw that a single chain undergoes a transition from a Luttinger liquid to 
The one-dimensional t — V model can be solved exactly via the Bethe 
ansatzE3. It is equivalent to two decoupled XXZ spin chains and so will undergo a metal-isulator transition a.t V — 2t. 
We note that the ring exchange process responsible for antiferromagnetic interactions discussed above will be absent 
in a single chain. Vojta, Hiibsch, and NoackH recently studied the model (^ on a ladder using the density matrix 
renormalisation group. They find that the charge ordered phase is destroyed for V < 2.5t but claim that there will be 
a charge gap for aM^V/t. The model in infinite dimensions was studied by Pietig, Bulla, and Blawid using dynamical 
mean-field theory.cll At low temperatures they found that for U = 2t, charge ordering occured for V > O.bt. 



III. ANTIFERROMAGNETIC EXCHANGE IN THE CHARGE-ORDERED PHASE 

For large V/t there will be charge ordering and the ground state will be an insulator with a charge gap of magnitude 
3V. There will be two possible ground states with the checker-board covering of the lattice (Fig. ||). This defines a 
new square lattice rotated by 45 degrees with respect to the original square lattice. It should be stressed that these 
ground states are distinct from the commensurate charge density waves, associated with a lattice distortion, and seen 
in some organics. To zero-th order in t/V all possible spin states will be degenerate. We consider a single placquette 
(Fig. 13) containing two spins. To second order in t/V both the singlet and triplet states have their energy lowered 
by — 4P/3F. The degeneracy of the singlet and triplet states is only broken to fourth order in t/V. We show below 
that this results in an effective antiferromagnetic exchange interaction 

At^ 

which acts along the diagonals of the original square lattice. Thus the spin degrees of freedom are described by a 
spin-i Heisenberg model on a square lattice. The Hamiltonian is 
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where Si denotes a spin operator on site i, and the sum (ij) runs over pairs of next-nearest neighbor lattice sites in 
the original square lattice. 

We now calculate the singlet-triplet splitting from fourth-order perturbation theory. If |5'o > is an eigenstate of 
^0 = ^X]<ij"> "t^ii^j then a perturbation Hi = H — Hq, which has no effect to third order, shifts the energy by 
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where the intermediate states labelled by {m, n,p} ^ are not degenerate with >■ The following process involving 
exchange of electrons around a placquette will contribute to a shift in the energy of the ground state. For the triplet 
states it can be represented as 
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The first, second, and fourth matrix elements are t and the third is — t. The intermediate states have energy 3V, 
AV , and 31^, respectively. (Note one needs to take into account the interaction with the neighbours not shown in the 
above representation). There are eight distinct ways of doing this exchange: at the first step there is four choices and 
at the third step there are two choices. Thus, the expression (||) implies that the triplet is increased in energy by 

A similar process for the singlet is 
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Thus, for the singlet this fourth-order process brings one back to the singlet wave function with a sign change and 
so the energy shift is opposite to that of the triplet. Hence, we arrive at (ph for the difference in energy between the 
singlet and triplet. 

It should be pointed out that there are also fourth order processes of the form |0 >— * \n >^ |0 >— > \p >^ |0 > 
which will produce a decrease in the ground state energy. However, because they produce the same shift for the 
singlet and triplet states we neglect them here. Our value of J is consistent with a recent calculationc2l of the effective 
exchange interaction in the charge ordered phase of the extended Hubbard model on a ladder when that result is 
rescaled to allow for different excitation energies of the intermediate states. For the ladder, the energies of the 
intermediate states are all 2V. Thus the ladder exchange is larger than for the square lattice by a factor of 9/2. 

We now consider whether this possible explanation for the origin of antiferromagnetism in the 9 type matacials is 
quantitatively reasonable. Roughly, it predicts that the value of J will be some fraction of t, typical valueseHl from 
Table | are of the order of 500-1000 K for the materials, with insulating gij©«nd states. 6'-(CiTET-TTF)2Br is an 
insulator with a charge gap of 600 meV and a value of ic3 of about 60 meVO Since the charge gap is 3V for t 
this gives V = 200 meV. Using J = At'^/dV^ gives J ~ 4 K (check) which is in reasonable agreement with the observed 
value of 6 K. We do not make quantitative comparisons of equation (|^) because they are not so clearly in the regime 
t < V, required for its validity. Fopsxample, for 6'-(BEDT-TTF)2Cu2(CN)[N(CN)2]2 Hiickel calculations give t - 80 
meV and the measured charge gapE3 is about 200 meV. This is inconsistent with the fact that the charge gap would 
be 3^ if t < V. 



IV. SLAVE BOSON THEORY FOR THE SU(N) VERSION OF THE EXTENDED HUBBARD MODEL 

WITH U oo 

We consider the SU(N) generalization of the Hamiltonian (^ for which the spin index, a, runs from 1 to N, and 
consider as a small expansion parameter assuming that N is large. Slave boson fields are introduced to allow 
treatment of the no double occupancy constraint required by the U —>■ oo limit. The effective action for the slave 
boson fields can be expanded in pccttets of The mean- field solution corresponds to the Gutzwiller approximation 
and is exact in the iV — > oo limit .OCj ThiSj-apBroach has been used to study other strongly correlated modeLs-such 
as the KondOjfflodel for magnetic imip«ritiesc3c3 in metals, the Hubbard modelo, the Hubbard-Holstein modeled and 
the AndersonO and Kondo lattices.E®! It has also been used to analyze the phase diagram of two-dimensional t-J 
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mode]lli~i3, and extended Hubbard model at half-filling.ii It is convenient-tp-jdescribe the projected Hilbert space 
associated with the Hamiltonian (||), using the slave boson representationSocJ The electron creation operator is 
replaced by: c-^ — fj^h, where the spinless charged boson operator, bi, is introduced to keep track of the empty sites, 
and fj^ is a fermion operator carrying spin. In order to preserve the anticommutation relation for the electrons the 
new operators must satisfy the local constraint 

flu + bjh = N/2. (8) 

For A'' = 2 it reduces to the psndition that either an electron or a boson can occupy each lattice site at all times. 
Following Kotliar and Liuo, we write the partition function in the coherent state path integral representation: 

Z = [ Db^DbDpDfDXcxpi - f L{T)dT\ (9) 



where the Lagrangian at imaginary time r is given by 

L{t) = E /-(^- - /^)/- + - ^ E iflfj^b]h + h.c.) 

i ij 



N 



(10) 



P = l/{kBT) at temperature T, and we have used the fact that c\^Cic — flfia- K is a static Lagrange multiplier 
enforcing the constraint (||). A sum from 1 to N is assumed whenever a repeated a index appears in the equations. 
The 1/N factors are introduced so that the Lagrangian is proportional to N. Tij = if i and j are nearest neighbours, 
Tij — t' \i i and j are next-nearest neighbours sitting along one of the diagonals of the square lattice, and Tij = 
otherwise. Vij = if i and j are nearest neighbours and is zero otherwise, /i is the chemical potential, which is 
fixed to give the average number of electrons per site, n =< f^fia >■ The conservation of the charge, q = N/2, is a 
consequence of a local U{1) symmetry because under the local gauge transformation: hi bie^^^^'^\ fia /io-e^^'^'^^ 
and Xi ^ Xi — dr6i{T), L{t) remains the same. To avoid possible infrared divergences it is convenient to choose the 
radial gauge where the boson amplitude becomes a real number, = and becomes a dynamical field: Ai(r). 
We introduce these fields in expression (|l^, and we use the relation flfia = N/2 — blbi to replace one pair of fermion 
operators in the V term so that we are left with a quadratic Lagrangian in the fermionic Grassmann variables. After 
integrating out the fermions, the effective Lagrangian for the boson fields becomes 

Lir) = ^{r,:(r)(a. + *A(r))r,(r) - *A.(r)^} 

i 

- NTr\n[{dr-fi + tX,{T) + 1 ^ - r,(r)n(r))H, - lr,;(r)r,,r,(r)] (11) 

I 



A. Mean-field solution 



The mean-field solution of the model is obtained by assuming that the boson fields are spatially homogeneous and 
time-independent: ri(r) = b and iXi{T) = X. The resulting free energy {F = —kBTln Z) is 

F''^{b, >^)^-jJ2 - lU^n) + \{b^ - y ) (12) 

where cj„ is a fermion Matsubara frequency. The mean-field eigenenergies are 

ek = ^Tk + A-M + 4y^ (13) 

with Tk = 2(cos(A:2;) + cos{ky) + ^ cos{kx + ky)) being the Fourier transform of Tij in units of the nearest-neighbour 
hopping t. 
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Minimization of the free energy (12) with respect to b and A gives 



A = ^/(ek)(trk + 4y). (14) 

k 

II is adjusted to give the correct electron fiUing, n = NJ2k /(^k) where /(e) is the Fermi-Dirac distribution function. 

The mean- field solution describes a renormalized Fermi liquid. The renormalization of the band is controlled by 6^, 
and the band is shifted from its bare position by A. The overall effect of the nearest-neighbours Coulomb interaction, 
V, at the mean-field level, reduces to a constant shift in the chemical potential. In the case of a quarter-filled band 
(n = 1/2) the bandwidth is reduced to half its bare value. The effective mass measured in magnetic oscillation 
experiments will then be increased by a factor of m* /m — — 2. Note that this is much smaller than the effective 
mass enham:ements that occur in materials described by the Hubbard model on the anisotropic triangular lattice at 
half-fiUing.El Therefore, for iV oo the quarter-filled t — V model behaves as a Fermi liquid with effective masses that 
are twice the bare ones. In the next subsection we consider the effect of the leading 1/N corrections. 



(16) 



B. Fluctuations about the mean- field solution 

We now consider how as V/t is increased, the Fermi liquid phase resulting from the mean-field solution becomes 
unstable to chareci ordering. The analysis is similar to the treatment of instabilities in the doped Hubbard model 
by Tandon et aZ.L3. We write the boson fields in terms of the static mean-field solution, (6, A), and the dynamic 
fluctuating parts: rj(r) ~ b + bSri{T), and iAi(r) = A -I- z<5Ai(r). Physically, Sri{T) is related to local fluctuations in 
the charge density. We substitute these expressions in (|lT|), introducing the Fourier transforms of Sri{T) and Xiir), 
and, expanding to second order in the boson fields, we obtain the effective action 

S = F'^^ + 5(2) (15) 

where the part of the action due to fluctuations in the boson fields is 

where Vn is a boson Matsubara frequency. The elements of the r(q, Vn) matrix are 

r.rA(q, l^n) = rAr(q, l^n) = N[— + l ^ ^ : (— r^(7k + Tfc+q) ^^q)] 

iv ek+q - Ek — iV 

TaaIq, i^n = -A^ > : — 17 

^ Ck+q - Ek - ll^n 

where Vk — '2{cos{kx) + cos{ky)) is the Fourier transform of Vij. Note that T\\ is the Lindhard function describing 
density-density fluctuations in the renormalised band structure. Since b^ is of order N (compare equation (|l^)), the 
above expressions show explicitly how the propagators of the boson fields -D(q, i/n) = r^^(q, t'n) are of order 0(1/A^), 
as they should. 



C. Charge ordering instability 

The condition for the stability of the Fermi liquid phase is that the quadratic form ( p^ is always positive. Then 
fluctuations in the charge density will increase the free energy. Since Faa > this is ensured if detr(q, ;/) > for 
all wavevectors q and frequency v. In order to find the critical value of V/t, which we shall denote (V/t)c, at which 
the system becomes unstable towards static charge ordering, we wish to find a q for which at some value of V/t, 
dctr = r,r.rr>;^ — F^rFrA = 0, at = 0. Tliis Condition reduces to: 
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[ ^ 2^ /(ek)(Tk-q - Tk + iV/t),Vk) + — (1 - V^) ] }^ 

k k 

2bH /(ek+q) - /(£k) ^rj. ,rp^,^b^_f. 

where A and &^ are the solution to the mean-field equations. 

We now concentrate on the case q —(t^, ti"), which is relevant to the formation of a charge ordered state at quarter 
filling in the 0-type materials (see Fig. For the square lattice case (t' — 0), fk+q = ^ Ck, and Eqn. ( |l8| ) reduces to 

where P(j(e) is the density of states at the Fermi surface per spin channel of the renormalized metal and W is the bare 
bandwidth of the metal. 

Before solving this equation numerically some insight can be gained by considering the case of a constant density of 
states. Taking a density of states of the form yOcr(e) = jr^jj, if —b'^W/2 < e < b'^W/2 and otherwise, eq.([l9|) can be 
simplified further. For this case, the critical ratio (V/t)c, at which (tt, tt) charge ordering occurs for a given electron 
band filling n is 

-^^y")^ -2n(f -l)ln(l-2f) 



«) + ite)_l]ln(l-2li) 

For N — 2, {V/t)c diverges at n = 0, n = 1 and n = 0.899. While the divergence at n = appears because it is 
not possible to have charge ordering when there is no charge in the lattice, the divergence at n = 1 is a consequence 
of the condition that there can be, at most, one electron at each lattice site: a charge ordered state of alternating 
singly and doubly occupied sites would cost infinite energy to be formed,. The divergence at n = 0.899 is non-trivial 
and presumably is a consequence of the finding made by Tandon et al£3 that close to half-filling, n > 0.88, the 1/N 
fluctuations drive the U oo Hubbard model into phase separation. Hence, the creation of a charge ordered state 
is forbidden by the breakdown of periodicity in the charge distribution of the lattice. At a quarter filling (n = 
equation ( po| ) gives {V/t)c = 0.78. This can be compared with the value of {V/t)c = 0.69 obtained from solving 
with the actual density of states for the square lattice. 



D. Effect of the diagonal hopping 



We can include the effect of a next-nearest neighbour hopping interaction in the analysis presented above. We have 
solved equation (|l8|), for different values of the t' /t ratio. Figure ^ shows how the critical value, {V/t)c, increases 
as the ratio, t' /t is increased. The critical value changes from 0.68 for the square lattice to {V/t)c ~ 0.95 for the 
triangular lattice {t' /t = 1). This can be understood if we plot the Fermi surfaces for different ratios of the hopping 
integrals for fixed band filling: n — 1/2 (see Fig. While for the square lattice case, t'/t — 0, the Fermi surface 
shows some remnants of the perfect nesting property present at half- filling, it gradually elongates along the ky = —k^ 
direction becoming elliptical as t'/t is increased£3 This effect makes it even harder to connect two points at the Fermi 
surface by the wavevector q = (tt, tt). The phase diagram in Fig. ^, also shows that it is possible to go from the 
metallic phase to the charge ordered state by varying either the t'/t or V/t ratios. Because the dependence of {V/t)c 
on t'/t is weak we conclude that the ratio V/t plays a more important role than t'/t in driving the metal-insulator 
transition, within the large-N approach used here. Hence, as a first approximation we are justified in neglecting the 
effect a£-the diagonal hopping, as was done in Section 
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MoriE3 found that the Coulomb repulsion V varied little with the angle 0. This is because it scales roughly with 
the inverse of the distance between the molecules. In contrast the hopping integrals tp and tc depend strongly on 9. 
tp varies by about a factor of five as 9 varies from 100 to 140 degrees. This is because the overlap integral depends 
exponentially on the distance between the molecules. Hence, the main effect of varying 9 is to change the band width. 
This is what will be driving the metal-insulator transition. j— , 

For 6I-(BEDT-TTF)2Cu2(CN)[N(CN)2]2 the measured charge gapEZi is about 200 meV, suggesting that V is of the 
order of 100 meV. Assuming that V does not vary much between materials the above calculations suggest that the 
critical value of the hopping integral t is about 100 meV. This is consistentc3 with the values in Table |; i.e., it is quite 
possible that the materials listed there are close to the metal-insulator transition, as is observed experimentally. 
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V. RELEVANCE TO /3"-(BEDT-TTF)2SF5YS03 



The family /5"-(BEDT-TTF)2SF5YS03 has been studicd0 with Y = CH2CF2, CH2, CHF. The first material is 
superconducting with a transition temperature of 5.2 K. Y = CH2 is insulating with a charge gap of 56 meV and 
evidence for charge ordering is found in the fact that altmiate molecules have a different bond length and phonon 
frequency associated with the central carbon double bondH The charges are estimated to be +0.6e and +0.4e where 
e is the electronic charge. Below room temperature the spin susceptibility decreases monotonically, consistent with 
a spin gap of 8 meV. Y=CHF is a bad metal and may undergo a metal-semiconductor transition below 10 K. It 
is estimated that alternate molecules have charges of +0A7e and -|-0.53e. A recent experimental studyES estimated 
charges of -|-0.43e and -|-0.57e in Y = CH2CF2. The Fermi surface of the metallic utiase of Y = CH2CF2 has been 
mapped out using angular-dependent magnetoresistance and magnetic oscillationsEj However, the meiallic phase 
differs significantly from a conventional Fermi liquid. First, in contrast tojiiost BEDT-TTF metalsjs even at a 
temperature as low as 14 K no Drude peak is present in optical conductivity.Ej Second, anomalous properties of the. 
magnetoresistance were recently interpreted in terms of a magnetic field induced superconductor-insulator transition.til 
The temperature dependence of the penetration depth in the superconducting phase was recently found to go like 
at low temperatures.E3 This is inconsistent with an s-wave state, but also deviates significantly from the linear 
temperature dependence expected for a d-wave state. 

The arrangement of the BEDT-TTF molecules within a layer of the family /3"-(BEDT-TTF)2SF5YS03 are shown 
in Figure Table |l| lists values of the hopping integrals calculated in the Hiickel approximation. Note that generally 
the diagonal terms a and a' are smaller than the vertical and horizontal terms. Hence, as a first approximation the 
system can be described by an anisotropic square lattice. However, we note that the main difference between the 
hopping parameters for the three different anions in Table || is that the diagonal terms a and a' vary as the anion is 
changed. This will change the proximity to the charge ordering instability and so lead to the three different ground 
states for these materials. 



VI. CONCLUSIONS 

In summary, we have argued that the essential physics of the electronic and magnetic properties of layered molecular 
crystals of the 9 type can be captured by an extended Hubbard model on the square lattice and at quarter filling. 
For large Coulomb repulsion {V » t) between nearest neighbours, the ground state is a charge-ordered insulator. 
Antiferromagnetic interactions arise due a novel fourth-order exchange process. A slave boson treatment was given 
of the SU(N) generalization of the model. It was found that for sufficiently small V/t the metallic phase was stable 
against charge ordering. _ 

We briefly mention the relation of this work to a recent paper of Mazumdar, Clay, and CampbelH. They have 
studied the extended Hubbard model at quarter filling on an anisotropic square lattice and discuss its relevance 
to a wide range of organics, but not those considered here. They argue that in the real materials the nearest- 
neighbour Coulomb repulsion is smaller than the critical value necessary to form the charge ordered state considered 
here. Coupling to phonons produces an insulating phase with a different kind of charge order (a bond-order- wave) . 
X-ray scattering experiments which can resolve the charge on individual molecules (due to different bond lengths) 
should be able to distinguish these two different ground states. There is, controversy about whether nucleaij-tnagnetic 
resonance measurements can distinguish these two charge orderings.ES The charge distribution observedtll for f3"- 
(BEDT-TTF)2SF5CH2S03 is consistent with charge ordering considered here. 

We acknowledge that the actual 9 type materials are more complicated than the simplest Hubbard model considered 
here. For example, along the diagonals pL the square lattice (corresponding to the vertical direction in Figure |l|) there 
is also Coulomlp-jepulsion. In fact, MoriE^ finds the corresponding V to be larger than along the horizontal and vertical 
directions. Seolia has shown how the latter can lead to competition between different charge ordered states (i.e., those 
associated with wave vector {tt,tt) and (0,7r)). Also, x-ray scattering suggests that in (BEDT-TTF)2RbM(SCN)4 
[M=Co,Zn] there is a structural transition associated with the charge ordering and that this changes the electronic 
structure in the insulating phase.E^I However, our view is that the t — V model on the square lattice captures the 
essential physics and first we need to understand it. 

Three outstanding questions concerning the t — V model at quarter filling need to be answered. 

(i) Is there superconductivity in the model? The idea that proximity to a quantum critical point increases— the 
tendency towards superconducting instabilities is supported by recent experiments on heavy fermion materials.o It 
was first shown by Scalapino, Loh, and HirschE3 that proximity to a spin-density wave or charge-density wave transition 
can lead to d-wave superconductivity. In a future publication-we will investigate whether charge fluctuations near the 
charge ordering transition can produce superconductivity.LjH 
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(ii) Are charge ordering, the charge gap, and antiferromagnetism destroyed at the same critical value of V jtl In 
Section III it was shown that for large V jt the ground state has a charge gap, charge ordering, and antiferromagnetism. 
In Section IV it was shown that, for V jt less than a critical value, the metallic phase is stable, at least in the large 
N limit. It is quite possible that the above three properties disappear at different values of Vjt. For the case 
of the quarter-filled extended Hubbard model on a ladder, numerical calculations found that the charge ordering 
disappeared below a non-zero value of Vjt, but the charge gap did not.E3 This unusual result may be an artefact 
of the one-dimensionality of the ladder. For the square lattice, this issue will probably be only resolved by careful 
numerical work. 

(iii) Does non-Fermi liquid behavior occur in the metallic phase near the quantum critical point? This is generally 
expectedEa and is observed in heavy fermion materialsE3 Slave boson theory has been used to show that in the doped 
Hubbard rpedel near a charge instability the quasi-particle scattering becomes singular leading to anomalous metallic 
properties O 
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FIG. 1. Arrangement of the donor molecules D (for example BEDT-TTF) within a layer for the &-D2X molecular crystals. 
The dashed rectangle denotes the unit cell. Typical values of the hopping integrals tp and tc are given in Table ^. Note that 
this geometry defines a tight binding model on an anisotropic triangular lattice which also can be viewed as a square lattice 
with hopping along one of the diagonals. The angle 6, and consequently the value of tp, varies with pressure or change in anion 
X. 
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FIG. 2. Schematic phase diagram of the 6'-(BEDT-TTF)2X and 6'-(BETS)2X famihesB0 showing competition between 
metallic, superconducting (SC), and charged ordered (CO) insulating phases. The horizontal axis is proportional to the angle 
(see Fig. 0) which is related to the hopping integral tp. Generally, increasing 6 decreases the bandwidth and so increases 
the importance of the electronic correlations. The vertical arrows denote the location of various materials at ambient pressure. 
The effect of pressure is to drive each material towards the right. 
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An antiferromagnetic interaction J occurs between spins along the diagonals. The spin degrees of freedom are described by the 
antiferromagnetic Heisenberg model on the square lattice. 
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FIG. 4. Phase diagram for the SU(N) version of the t ~ t' — V model at quarter filling and zero temperature, to leading 
order in 1/N with N=2. This shows that the diagonal hopping t' {tc in Figure hi) has little effect on the critical value of V/t at 
which the metallic phase becomes unstable to the charge ordering shown in Figure ^. 
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crystals. Typical values of the hopping integrals are given in Table In some of the materials the unit cell is larger and so 
the hopping integrals can have two values. 
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TABLE I. Hopping integrals for various 9-type crystals calculated by the Hiickel method. Two values are given for the case 
where the unit catt is larger. The temperature of the ipffital-insulator transition, Tmi, is also given. Note the general trend, 
observed by MoricJ, that as tp increases Tmi decreases.E3 



Material 


tp (moV) 


tc (mcV) 


Reference 


Tmi (K) 


(BEDT-TTF)2l3 


42 


64 




30 






(BETS)2Ag(CN)2 


392,398 


-1,38 








< 4 


(BETS)4Cu2Cl6 


380-467 


-12 - -1-56 




73 




< 4 


(BETS)2CsCo(SCN)4 


366 


-2 




12 




10? 


(BETS)2CsZn(SCN)4 


372 


-10 




12 




10? 


(BETS)2RbCo(SCN)4 


382 


-72 




12 




20 


(BETS)2RbZn(SCN)4 


347 


-46 








? 


(BEDT-TTF)2CsCo(SCN)4 


106 


-5 




29 




20 


(BEDT-TTF)2CsZn(SCN)4 


108 


-10 




29 




20 


(BETS)4TaF6 


-30 


110 




74 




70 


(BEDT-TTF)2RbCo(SCN)4 


99 


-33 




29 




190 


(BEDT-TTF)2RbZn(SCN)4 


94 


-24 




29 




190 


(BEDT-TTF)2Cu2(CN)[N(CN)2]2 


79 


-30 




27 




220 


(BEDT-TTF)2TlCo(SCN)4 


100 


-48 




29 




250 


(BDT-TTP)2Cu(NCS)2 


-86,-91 


-41 




17 




250 


(CiTET-TTF)2Br 


-54,-43 


-58 




17 




> 300 



TABLE IL Hopping integrals in Figure || for the family /3"-(BEDT-TTF)2SF5YS03 calculated by the extended Hiickel 
method in Reference uA. At low temperatures the materials are a superconductor, a bad metal, and a charge ordered insulator, 
respectively. 



Y 


c (meV) 


d (meV) 


a (meV) 


a' (meV) 


CH2CF2 


260 


140 


120 


55 


CHF 


260 


130 


35, 86 


95, 100 


CH2 


260 


120 


85 


12 



15 



